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Abstract 

A natural way to deal with multiple, partially conflicting objectives is turning all the objectives 
but one into budget constraints. Some classical polynomial-time optimization problems, such 
as spanning tree and forest, shortest path, (perfect) matching, independent set (basis) in a ma- 
troid or in the intersection of two matroids, become NP-hard even with one budget constraint. 
Still, for most of these problems deterministic and randomized polynomial-time approximation 
schemes are known. In the case of two or more budgets, typically only multi-criteria approx- 
imation schemes are available, which return slightly infeasible solutions. Not much is known 
however for the case of strict budget constraints: filling this gap is the main goal of this paper. 

We show that shortest path, perfect matching, and spanning tree (and hence matroid basis 
and matroid intersection basis) are inapproximable already with two budget constraints. For 
the remaining problems, whose set of solutions forms an independence system, we present de- 
terministic and randomized polynomial-time approximation schemes for a constant number k 
of budget constraints. Our results are based on a variety of techniques: 

1. We present a simple and powerful mechanism to transform multi-criteria approximation 
schemes into pure approximation schemes. This gives, for example, deterministic approx- 
imation schemes for fc-budgeted forest and bipartite matching, and randomized approxi- 
mation schemes for /c-budgeted matching, independent set in matroids, and independent 
set in the intersection of two representable matroids. 

2. We show that points in low dimensional faces of any matroid polytope are almost integral, 
an interesting result on its own. This gives a deterministic approximation scheme for 
fc-budgeted matroid independent set. 

3. We present a deterministic approximation scheme for 2-budgeted matching. The backbone 
of this result is a purely topological property of curves in R 2 . 



1 Introduction 

In many applications, one has to compromise between several, partially conflicting goals. Multi- 
Objective Optimization is a broad area of study in Operations Research, Economics and Computer 
Science [bH 129] , A variety of approaches have been employed to formulate such problems 
including Goal Programming [3], Pareto-Optimality [10], and Multi-Budgeted Optimization [29J. 
We adopt the latter approach and cast one of the goals as the objective function, and the others as 
budget constraints. More precisely, we are given a (finite) set T of solutions for the problem, where 
each solution is a subset S of elements from a given universe E (e.g., the edges of a graph). We 
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are also given a weight function w : T — > Q + and a set of k = 0(1]_| length functions l{ : T —> Q + , 
1 < i < k, that assign a weight w(S) := ^eeS w ( e ) an< ^ an zth-length £i(S) := ^ ee s^i(e), 
1 < i < k, to every candidate solution S. For each length function £i, there is a budget Lj € Q + . 
The multi-budgeted optimization problem can then be formulated as follows: 

maximize/minimize w(S) subject to S £ J-, £i(S) < Lj, 1 < i < k. (1) 

We next use OPT to denote an optimum solution. 

Following the literature on the topic, we focused on the set of problems below: 

• A;-budgeted (perfect) matching: T is given by the (perfect) matchings of an undirected 
graph G= (V,E). 

• /c-budgeted SPANNING tree (forest): T is given by the spanning trees (forests) of G. 

• /c-budgeted shortest path: T is given by the paths connecting two given nodes s and t 
in G. 

• /c-budgeted matroid independent set (basis): T is given by the independent sets 
(bases) of a matroid M = (E,lj^. 

• A;-budgeted matroid intersection independent set (basis): T is given by the inde- 
pendent sets (bases) in the intersection of two matroids M\ = (E,Ii) and Mi = (E,l2). 

All the problems above are polynomial-time solvable (see, e.g., [31 ) in their unbudgeted version 
(k = 0), but become NP-hard [U El [12] even for a single budget constraint (k = 1). For the case 
of one budget (k = 1), polynomial-time approximation schemes (PTASs) are known for spanning 
tree [28] (see also [T7j), shortest path [32] (see also [13 [21]), and matching [6] (see also [5]). 
The approach in [28J easily generalizes to the case of MATROID BASIS. A PTAS is also known for 
MATROID intersection independent set [6j. The results in [6] do not generalize to the case of 
PERFECT MATCHING and MATROID INTERSECTION BASIS. No approximation algorithm is known 
for the problems above in the case k > 2 (excluding multi-criteria algorithms which provide slightly 
infeasible solutions): investigating the existence of such algorithms is the main goal of this paper. 

1.1 Our Results 

We start by observing that several of the mentioned problems are inapproximable already for two 
budget constraints. More precisely, the corresponding feasibility problem is NP-complete. The 
simple proof of the following theorem is given in the appendix. 

Theorem 1. For k > 2, unless P = NP there is no polynomial-time approximation algorithm for 
fc-BUDGETED SPANNING TREE, fe-BUDGETED SHORTEST PATH, fe-BUDGETED PERFECT MATCHING, 
fc-BUDGETED MATROID BASIS, and fc-BUDGETED MATROID INTERSECTION BASIS. 

The remaining problems have a common aspect: the set of solutions J- forms an independence 
system. In other terms, for S € T and S' C S, we have S' £ T . For these problems, we present 
deterministic and randomized approximation schemes, based on a variety of techniques. 

x The assumption that k is a constant is crucial in this paper. 

2 We recall that E is a finite ground set and I C 2 E is a nonempty family of subsets of E (independent sets) 
which have to satisfy the following two conditions: (i) I £ X, J C I => J £ I and (ii) I, J £ I, |/| > \J\ => 3z G 
I \ J : J U {z} G I. A basis is a maximal independent set. For all matroids used in this paper we make the usual 
assumptions that independence of a set can be checked in polynomial time. For additional information on matroids, 
see e.g. [5T| . 
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Our first result (see Section [2]) is a simple but powerful mechanism to transform a multi-criteria 
PTAS, i.e. a PTAS that might violate the budgets by a small multiplicative factor, into a pure 
PTAS, where no budget is violated. Similarly, a multi-criteria polynomial randomized-time ap- 
proximation scheme (PRAS) can be transformed into a pure PRAS. 

Theorem 2. (Feasibilization) Let Vi n d be a k-budgeted problem where the set of solutions T 
is an independence system. Suppose that we are given an algorithm A which, for any constant 
5 > 0, computes in polynomial time an (1 — 6) (resp., expected (1 — 6)) approximate solution to 
Find violating each budget by a factor at most (1 + 5). Then there is a PTAS (resp., PRAS) for 

Find- 

The idea behind the proof is showing that a good solution exists even if we scale down the budgets 
by a small factor. This is done by applying a greedy discarding strategy similar to the greedy 
algorithm for knapsack. Applying a multi-criteria PTAS (given as a black box!) to the scaled 
problem gives a feasible solution for the original one, of weight close to the optimal weight. 
To the best of our knowledge, this simple result was never observed before. Indeed, it implies im- 
proved approximation algorithms for a number of problems. A general construction by Papadim- 
itriou and Yannakakis |25j provides multi-criteria PTASs (resp., PRASs) for problems whose exact 
version admits a pseudo-polynomial-time (PPT) deterministic (resp., Monte-Carlo) algorithm. We 
recall that the exact version of a given optimization problem asks for a feasible solution of exactly 
a given target weight. Combining their approach with our mechanism one obtains approximation 
schemes for several problems. For example, using the PPT-algorithm for exact forest in [3], one 
obtains a PTAS for /c-budgeted forest. Similarly, the Monte-Carlo PPT-algorithm for exact 
matching in [21] gives a PRAS for /c-budgeted matching. The Monte-Carlo PPT-algorithms 
for exact MATROID intersection independent set in [8], which works in the special case of 
representable matroida^l, implies a PRAS for the corresponding budgeted problem. 
Of course, one can also exploit multi-criteria approximation schemes obtained with different tech- 
niques. For example, exploiting the multi-criteria PTAS in [13] for /c-budgeted matching in 
bipartite graphs, which is based on iterative rounding, one obtains a PTAS for the same prob- 
lem. Very recently [9], a multi-criteria PRAS for /c-budgeted MATROID independent set, based 
on dependent randomized rounding, has been presented. This implies a PRAS for /c-budgeted 
MATROID INDEPENDENT SET. 

Corollary 3. There are PTASs for /c-budgeted forest and /c-budgeted matching in bipartite 
graphs. There are PRASs for /c-budgeted matching, /c-budgeted matroid independent set, 
and Zc-budgeted matroid intersection in representable matroids. 

Based on a different, more direct approach, we are able to turn the PRAS for /c-budgeted matroid 
independent set into a PTAS. The main insight is the following structural property of faces of 
the matroid polytope which might be of independent interest (proof in Section [3]). 

Theorem 4. Let M = (E,I) be a matroid and let F be a face of dimension d of the matroid 
polytop^ Px- Then any x € F has at most 2d non-integral components. Furthermore, the sum of 
all fractional components of x is at most d. 

A PTAS can then be easily derived as follows. We first guess the k/e elements Eh of largest weight 
in the optimum solution in a preliminary phase, and reduce the problem consequently. This guessing 

3 A matroid M = (S, I) is representable if its ground set S can be mapped in a bijective way to the columns of a 
matrix over some field, and I C S is independent in M iff the corresponding columns are linearly independent. 

4 For some given matroid M = (E,X), the corresponding matroid polytope Px is the convex hull of the incidence 
vectors of all independent sets. 
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step guarantees that the maximum weight w max of an element in the reduced problem satisfies 
^w max < ew(Eh)- For the reduced problem, we compute an optimal fractional vertex solution x* 
to the LP which seeks to find a maximum weight point in the matroid polytope intersected with the 
k budget constraints. Since x* is chosen to be a vertex solution, and only k linear constraints are 
added to the matroid polytope, x* lies on a face of the matroid polytope of dimension at most k. We 
then round down the fractional components of x* to obtain an incidence vector x which corresponds 
to some independent set El- By Theorem HI \x* — x\ < k, and hence, w(El) > w(x*) — /cw max . 
Then, it is not hard to see that Eh U El is a (1 — e)-approximate feasible solution for the starting 
problem. 

Corollary 5. There is a PTAS for fc-BUDGETED matroid independent set. 

Eventually, we present a PTAS (rather than a PRAS as in Corollary [3]) for 2-budgeted matching 
(see Section 

Theorem 6. There is a PTAS for 2-budgeted matching. 

Our PTAS, works as follows. Let us confuse a matching M with the associated incidence vector xm- 
We initially compute an optimal fractional matching x*, and express it as the convex combination 
x* = a.\X\ + CK2X2 + a 3 x 3 °f three matchings x%, X2, and X3. Then we exploit a patching procedure 
which, given two matchings x' and x" with high Lagrangian weight and a parameter /i € [0,1], 
computes a matching z which is not longer than x^ := fix' + (1 — /j,)x" with respect to both lengths, 
and has a comparable weight. This procedure is applied twice: first on the matchings x\ and X2 
with parameter [i = a\j{a\ + 02), hence getting a matching z' . Second, on the two matchings z 1 
and X3 with parameter /i = («i + 02) /(ai + «2 + 03). The resulting matching z" is feasible and 
almost optimal (modulo a preliminary guessing step). 

Our patching procedure relies on a topological property of curves in M 2 , that we prove via Jordan's 
curve theorem [22]. An extension of the property above to curves in M fc would imply a PTAS for 
/c-budgeted matching: this is left as an interesting open problem. 

1.2 Related Work 

There are a few general tools for designing approximation algorithms for budgeted problems. 
One basic approach is combining dynamic programming (which solves the problem for polyno- 
mial weights and lengths) with rounding and scaling techniques (to reduce the problem to the 
case of polynomial quantities). This leads for example to the FPTAS for 1-budgeted shortest 
path [161 [21] [32] . Another fundamental technique is the Lagrangian relaxation method. The basic 
idea is relaxing the budget constraints, and lifting them into the objective function, where they are 
weighted by Lagrangian multipliers. Solving the relaxed problem, one obtains two or more solutions 
with optimal Lagrangian weight, which can - if needed - be patched together to get a good solution 
for the original problem. Demonstrating this method, Goemans and Ravi [28] gave a PTAS for 1- 
budgeted SPANNING tree, which also extends to 1-BUDGETED MATROID BASIS. Using the same 
approach, with an involved patching step, Berger, Bonifaci, Grandoni, and Schafer [6] obtained a 

PTAS for 1-BUDGETED MATCHING and 1-BUDGETED MATROID INTERSECTION INDEPENDENT SET. 

Their approach does not seem to generalize to the case of multiple budget constraints. 
The techniques above apply to the case of one budget. Not much is known for problems with two 
or more budgets. However, often multi- criteria approximation schemes are known, which provide 
a (1 — e)-approximate solution violating the budgets by a factor (1 + e). First of all, there is 
a very general technique by Papadimitriou and Yannakakis [25j, based on the construction of e- 
approximate Pareto curves. Given an optimization problem with multiple objectives, the Pareto 
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curve consists of the set of solutions S such that there is no solution S' which is strictly better 
than S (in a vectorial sense). Papadimitriou and Yannakakis show that, for any constant e > 0, 
there always exists a polynomial-size e-approximate Pareto curve A, i.e., a set of solutions such 
that every solution in the Pareto curve is within a factor of (1 + e) from some solution in A on 
each objective. Furthermore, this approximate curve can be constructed in polynomial time in the 
size of the input and 1 /e whenever there exists a PPT algorithm for the associated exact problem. 
This implies multi-criteria FPTASs for /c-budgeted spanning tree and /c-budgeted shortest 
path. Furthermore, it implies a multi-criteria FPRAS for /c-budgeted (perfect) matching. 
The latter result exploits the Monte-Carlo PPT algorithm for exact matching in [23]. Our PRAS 
improves on these results, approximation-wise (the running time is larger in our case). 
Recently, Grandoni, Ravi and Singh |13] showed that the iterative rounding technique is an al- 
ternative way to achieve similar (or better) results. The idea behind iterative rounding |18] (see 
also, e.g., [2J 20J) is to consider a linear relaxation of the problem, compute an optimal fractional 
solution, and round one of its variables. The process is then iterated on the residual problem until 
a feasible integral solution is obtained. This approach can be enhanced with a relaxation step, 
where a constraint which cannot be violated too much is relaxed (i.e., deleted). Using this method, 
Grandoni et al. obtain a multi-criteria PTAS for /c-budgeted SPANNING tree, which computes a 
solution of optimal cost violating each budget by a factor (1 + e). This improves, approximation- 
wise, on the result in [25] for the same problem (where the solution returned is suboptimal). The 
authors also show how to obtain a deterministic (rather than randomized [25j) multi-criteria PTAS 
for A;-BUDGETED MATCHING in bipartite graphs. 

All the mentioned problems are easy in the unbudgeted version. Given an NP-hard unbudgeted 
problem which admits a p approximation, the parametric search technique in [23J provides a multi- 
criteria kp approximation algorithm violating each budget by a factor kp for the corresponding 
problem with k budgets. Other techniques lead to logarithmic approximation factors (see, e.g., 

P3I251E7J). 

2 A Feasibilization Mechanism 

In this section we illustrate our feasibilization mechanism, proving Theorem [2J 

Proof of Theorem^ Let e £ (0, 1] be a given constant. Consider the following algorithm. Initially 
we guess the h = k/e elemental Eh of OPT of largest weight, and reduce the problem conse- 
quently, hence getting a problem V' . Then we scale down all the budgets by a factor (1 — 5), 
and solve the resulting problem V" by means of A, where 5 = s/(k + 1). Let El be the solution 
returned by A. We eventually output Eh U El- 

Let OPT' and OPT" be the optimum solution to problems V' and V" , respectively. We also denote 
by \J i and L" the ith budget in the two problems, respectively. Eventually, let w max be the largest 
weight in V' and V" . We observe that trivially 

(a) w(OPT) = w(E H ) + w(OPT') and (b) w max < w(E H )/h. 

Let us show that 

w(OPT") > w(OPT')(l - k5) - kw max . (2) 
5 To avoid inessential technicalities, we assume 1/e 6 N. 

6 As usual, by reducing we mean decreasing each budget Li by ti(Eu) and removing all the elements of weight 
strictly larger than mhiegE^ w(e). By guessing we mean trying all the 0(m h ) subsets of h elements. 
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Consider the following process: for each length function i, we remove from OPT' the element e 
with smallest ratio w(e)/£j(e) until ii(OPT') < (1 — 5)1^. Let E{ be the set of elements removed. 
It is not hard to see that w(-EA < 5w(OPT') + w max . It follows that OPT' — Uj-Ej is a feasible 
solution for V" of weight at least w(OPT')(l — 5k) — kw max , proving Q. 
We observe that El is feasible for V' since, for each i, 

£i(E L )<(l + 6)L'! = (l + S)(l-6)L' i <L' i . 

As a consequence, the returned solution Eh U El is feasible. Moreover, when A is deterministic, 
we have 

w(E H ) + w(E L ) > w(E H ) + (1 - 6)w(OPT") > w(E H ) + (1 - <5)(u;(OPT')(l - 5k) - kw max ) 

0) 

> (1 - k/h)w(E H ) + (1 - <5(/c + l))w(OPT') > (1 - e)(w(E H ) + w(OPT')) 
= {l-e)w{OPT). 

The same bound holds in expectation when A is randomized. □ 

3 A PTAS for /c-Budgeted Matroid Independent Set 

It is convenient to consider weights w and lengths ii as vectors in Q E . We denote by I the matrix 
whose ith column is £i, and let L = (Li, . . . ,L^) T . To every matroid M = (E,I), a rank function 
r : 2 E — >• N is associated; it is defined by r(S) = max{| J\ \ J C S, J € X}. The matroid polytope Pi 
is the convex hull of the characteristic vectors \i °f the independent sets I E T and is described by 
the following set of inequalities (see [31] for more details): 

P x = conv{x/ :lel} = {x>0: x{S) < r(S) VS C E}. 

Proof of Theorem^ Let m = \E\. We assume that the matroid polytope has full dimension, i.e., 
dim(Pj) = m, which is equivalent to saying that every element e € E is independent. This can 
be assumed without loss of generality since if {e} ^ I for some e £ E, then we can reduce the 
matroid by deleting element e. Since dim(Pj) = m and dim(i ? ) = d, F can be described by the 
inequality system of Px, where m — d linearly independent inequalities used in the description of 
Px are turned into equalities. More precisely, there are iV C E and A\, . . . , C E such that 

F = {x € P x | x(e) = Ve € A, zr(A) = r(Ai) V< G {1, . . . , A;}}, 

and |A| + = m — d. By standard uncrossing arguments, we can assume that the sets Ai form a 
chain, i.e., Ai C ^4 2 C • • • C A^ (see for example [15] H] for further information on combinatorial 
uncrossing). We prove the claim by induction on the number of elements of the matroid. The 
theorem clearly holds for matroids with a ground set of cardinality one. First assume N ^ and 
let e € N. Let M' be the matroid obtained from M by deleting e, and let F' be the projection of 
F onto the coordinates corresponding to N \ {e}. Since F' is a face of M' , the claim follows by 
induction. Henceforth, we assume N = which implies k = m — d. Let Aq = and Bi = Aj \ Aj_x 
for i 6 {1, . . . , k}. In the following we show that we can assume 

< r(Ai) - r(Ai_i) < V i G {1, . . . , fc}. (3) 

Notice that < r(Aj) — r(Aj_i) < \P>i\ clearly holds by standard properties of rank functions 
(see [3T| for more details). Assume that there is i G {!,...,&} with r(A«) = r(Aj_i). Since all 
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points x E F satisfy x(Ai) = r{Ai) and x(A;_i) = r(^4j_i), we have x(Bj) = 0. Hence for any 
e E -Bj, we have x(e) = for x E i* 1 . Again, we can delete e from the matroid, hence obtaining a 
smaller matroid for which the claim holds by the inductive hypothesis. Therefore, we can assume 
r(Ai) > r(Ai-i) which implies the left inequality in 0. 

For the right inequality assume that there is i E {1, . . . , k} with r{Ai) — r(Aj_i) = \B{\. Hence, 
every x E F satisfies x(Bi) = \Bi\, implying x(e) = 1 for all e £ ft. Let e E Bi, and let F' be 
the projection of the face F onto the components N \ {e}. Since F' is a face of the matroid M' 
obtained from M by contracting e, the result follows again by the inductive hypothesis. 
Henceforth, we assume that ([3]) holds. This implies in particular that \B{\ > 1 for i E {1, . . . , k}. 
Since Yli=i — m > we have k < m/2, which together with k = m — d implies d > m/2. The 
claim of the theorem that x E F has at most 2d non-integral components is thus trivial in this case. 
To prove the second part of the theorem we show that if ([3]) holds then x(E) < d for x E F. For 
igFwe have 



:{E) = x(E \A k ) + Y, x(Bi) < \E\ - \A k \ + ^(r(^) - r(^_0) 

i=l i=l 

< \E\ -\A k \ + ^2(\Ai\ - \Ai- X \ - 1) = m - k = d, 



where the first inequality follows from x(E \ A k ) < \E\A k \ and x(B{) = r{Ai) — r(^4j_i), and the 
second inequality follows from ([3|). □ 



4 A PTAS for 2-BUDGETED MATCHING 

In this section we present our PTAS for 2-budgeted matching. We denote by M the set of 
incidence vectors of matchings. With a slight abuse of terminology we call the elements in A4 
matchings. Let Pm be the matching polyhedron. Analogously to Section [3l let I = (^1,^2) and 
L = (Li,L2) r . A feasible solution in this framework is a matching x E M. such that £ T x < L. For 
two elements z',z" E [0,1]^, we define their symmetric difference z'Az" E [0, 1] E by (z' Az")(e) = 
\z'(e) — z"{e)\ for all e E E. In particular, if z' and z" are incidence vectors, then their symmetric 
difference as defined above corresponds indeed to the symmetric difference in the usual sense. Recall 
that, when z' and z" are matchings, z'Az" consists of a set of node-disjoint paths and cycles. 
We start by presenting a property of curves in R 2 (Section 14. ip . This property is used to derive 
the mentioned patching procedure (Section I4.2p . Eventually, we describe and analyze our PTAS 
(Section H3|). 

4.1 A Property of Curves in R 2 

We next describe a topological property of polygonal curves in M?, which will be crucial in our 
proofl A curve in M. 2 is a continuous function / : [0,r] — > M 2 for some r E M+. A curve is called 
polygonal if it is piecewise linear. 

For a E [0, r], let f a : [0, r] ->■ M 2 be the following curve. 

ra) = {f {t + a ^- + /(°) if 1 + a < r ' 

\f(r) ~ f(a) + f(a + t-r) if t + a > r. 

7 The lemma even holds for general (non-polygonal) curves. However, since we only need polygonal curves in our 
setting we restrict ourselves to this case since it simplifies the exposition. 
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Observe that / a (0) = /(0) and f a (r) = /(r) for any a G [0, r]. The next lemma shows that any 
point x on the segment between /(0) and /(t) is contained in some curve f a . 

Lemma 7. Let f : [0, t] — >• R 2 6e a polygonal curve, and let fi G [0, 1]. T/ien i/iere are a,t £ [0, r] 
siic/i i/wrf / a (t) = /i/(0) + (1 - //)/(r). 

We next give an intuitive description of the proof of the lemma: a formal proof is given in the 
appendix. Let / = (/l,/^)- Since the statement of the lemma is independent of changes in the 
coordinate system (and the claim is trivial for /(0) = /(r)), we can assume that /(0) = (0,0) and 
/(t) = (r, 0) for some r > 0. The Gasoline Lemma [6] states that there is a\ G [0, r] such that 
f^ft) — Vt G [0,t]. In particular, this condition is satisfied by choosing a\ G argmin{/2(t) | i € 
[0, r]}. Analogously, for 02 G argmax{/2(i) | t G [0,r]}, /^(i) < Vi G [0, t\. Hence, we have two 
curves, f ai and / a2 , one above and the other below the x-axis, both with the same endpoints (0, 0) 
and (r, 0). (See Figure [TJ. Furthermore, for a ranging from a\ to 02 (in a circular sense), the curve 
f a continuously transforms from f ai to f a2 , always maintaining the same endpoints. Then it is 
intuitively clear that the union of the curves f a spans all the points on the segment from (0, 0) to 
(r, 0), hence proving the claim. 



4.2 The Patching Procedure 

In this section we describe a patching procedure which, given two matchings x' and x" and a 
parameter [i G [0, 1], computes a matching z satisfying t T z < £ T x^, where x^ := \ix' + (1 — fj)x" is 
a convex combination of the first two matchings. Furthermore, the weight w T z is close to w T x ll , 
provided that x' and x" have a sufficiently large Lagrangian weight, which is defined as follows. Let 
A^, A2 £ be a pair of optimal dual multipliers for the budgets in the linear program max{w T x | 
x G Pm,£ T x < L}. The Langrangian weight of x G [0,1] E is C(x) = w T x — (A^A^X^x - L). 
Notice, that by the theory of Lagrangian duality we have w* = m&x{£(x) \ x G -Pm}> where w* is 
the weight of an optimal LP solution, i.e., w* = max{w T x | x G P_m,£ t x < L} (see [19! for more 
information on Lagrangian duality). 
We need the following notion of almost matching. 

Definition 8. For r G N, an r-almost matching in G is a (possibly fractional) vector y G [0, 1] E 
such that it is possible to set at most r components of y to zero to obtain a matching. 

We denote by A4 r the set of all r-almost matchings in G. Given an r-almost matching y, we let a 
corresponding matching z G Ai be a matching obtained by setting to zero the fractional components 
of y, and then computing a maximal matching in the resulting set of edges (in particular, we might 
need to set to some 1 entries of y to obtain z). Notice that w T z > w T y — rw max , where w m&x is 
the largest weight. 

Our patching procedure first constructs a 2-almost matching y, and then returns a corresponding 
matching z. We next show how to compute y. Let us restrict our attention to the following set of 
candidate 2-almost matchings. Recall that s = x'Ax" is a set of paths and cycles. We construct an 
auxiliary graph C, consisting of one cycle (eo, ei, . . . , e r _i), with the following property: there is a 
bijective mapping between the edges of C and the edges of s such that two consecutive edges of C 
are either consecutive in some path/cycle or belong to different paths/cycles. This can be easily 
achieved by cutting each cycle, appending the resulting set of paths one to the other, and gluing 
together the endpoints of the obtained path. For t G [0, r], we define s(t) G [0, 1] E as 



( S (t))(e) 



1 if e = ej, i < [t\ ; 

t— \ t\ if e = ej, i= [t\; 
otherwise. 
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Moreover, for a, t £ [0, r], we define 
[0,l] s 9^(i) = 



s(a + t)— s(a) if a + t < r; 

s(a + i — r) + s(t) — s(a) if a + t > r. 



Intuitively, a and (a + i) (mod r) define a (fractional) subpath of C, and s a (i) is the (fractional) 
incidence vector corresponding to that subpath. Eventually we define 

y a (t) :=x'As a (t). 

Note that y a (t) is equal to x' and x" for t = and t = t, respectively. 
Lemma 9. For any a,t G [0,r], y a (i) is a 2-almost matching. 

Proof. One can easily observe that a matching can be obtained by setting the two components of 
y a {t) to zero that correspond to the edges e^ a j and e [(«_)_*) (mod r)J- d 

The following lemma shows that, in polynomial time, one can find a 2-almost matching y with 
lengths £ T y equal to the lengths of any convex combination of the two matchings x' and x". 

Lemma 10. Let \i € [0,1] and x^ = fix 1 + (1 — fi)x" . In polynomial time, a,t € [0,r] can be 
determined such that £ T y a (t) = £ T x^. 

Proof. Let / : [0,r] — > M 2 be the polygonal curve defined by f(t) = £ T y°(t). Since /(0) = £ T x' 
and /(t) = £ T x", we have by Lemma[T]that there exists a,t 6 [0,r] such that / a (i) = £ T x^. Since 
/ a (i) = £ T y a (t), y := y a (i) satisfies the claim. 

The values of [a\ and [ a + ^J can be guessed in polynomial time by considering 0(n 2 ) possibilities. 
Given those two rounded values, the actual values of a and t can be obtained by solving a linear 
program with a constant number of variables and constraints. □ 

Our patching procedure simply computes a 2-almost matching y = y a {t) with £ T y = (Fxn, exploit- 
ing the lemma above, and then returns a corresponding matching z, by applying the procedure 
explained in the proof of Lemma [9l Trivially, £ T z < £ T y = £ T x^. We next show that, if x' and x" 
have sufficiently large Lagrangian weight, then the weight of z is close to the weight of x^. 

Lemma 11. Assume C(x') > w* — T and C{x") > w* — T for some T € R + . Then the matching z 



returned by the patching procedure satisfies w T z > w T x^ — 2w max — T and £ T z < £ T x,_ 



Proof. By Lemma [10] we have i T y = £ T x^, and since z < y, we get £ T z < £ T x^. Let x^ 



x' + x" — x^ = (1 — /j,)x' + fix". Since C{x') > w* — T, C{x") > w* — T and C is linear, we have 
> w* — r and C{x,j) > w* — V. Recall that y = y a (t) for a proper choice of a, t € [0, r]. Let 
y := x' + x" — y. Notice that y = y a (r — t) where a' = (a + i) (mod r), and hence, y is also a 
2-almost matching by Lemma GO Let ~z be the matching corresponding to y obtained by applying 
the procedure explained in the proof of Lemma [S] to y. Notice that the pairs (z, y) and (z, y) differ 
on the same two (or less) components. Hence 

rp rp rp rp rp rp 

w z + w z + 2w max >w y + w y = w x^ + w x^. (4) 

Since y + y = x^ + x M and £ T y = l T x^, we get £ T y = £ T x M . Thus, £ T z < £ T x^ since z <y. This 
can be rewritten as £(z) — w T z > Ci^x^) — w 7 ^. Since C(x^) > w* — T and C(z) < w*, we obtain 
w T z < w T x^ + r. Combining this result with (JH) implies w T z > w T x^ — 2w max — T. □ 
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4.3 The Algorithm 



Our PTAS works as follows. Initially it guesses the 6/e heaviest edges Eh in the optimum solution, 
and reduces the problem consequently. Then it computes a vertex x* G Pm of the polytope 
{x € Pm I £ T x < L} of maximum weight w* := w T x*. As x* is a vertex solution of the polytope 
Pm with two additional constraints, it lies on a face of Pm of dimension at most two. Hence, by 
Caratheodory's Theorem, x* can be expressed as a convex combination x* = a\X\ + 02^2 + 03^3 
of three matchings xi,x 2 ,x 3 G Pm- Let // = a\/(a\ + a 2 ) and fi" = (a\ + a 2 )/(a\ +02 + 03). 
Applying Lemma [TT1 to xi and X2 with /x = //, a matching z' is obtained. Applying Lemma [TT1 to 
2/ and x 3 with u = fi", we obtain a matching z". The algorithm returns z" plus Eh- 

Proof of Theorem Consider the algorithm above. The initial guessing can be performed in 
0{\E\ % / £ ) time. Since it is possible to efficiently separate over Pm, x* can be computed in poly- 
nomial time |31j. The same holds for the decomposition of x* into three matchings by standard 
techniques (see for example [3Qj ) - Lemma [TOl implies that the patching can be done in polynomial 
time. Hence the proposed algorithm runs in polynomial time as claimed. 

Since C(x*) = w* and C(x) < w* for x € Pm, we get C(x\) = C(x 2 ) = £(x 3 ) = w*. Let 
u := /i'xi + (1 — \J)x2 and v := [i"z' + (1 — /i")^3- By Lemma[TTl matching z' satisfies l T z' < £ T u 
and w T z' > w T u — 2w max . Since u is a convex combination of x\ and X2, we have C{u) = w*. 
Furthermore, by the relations between the lengths and weight of z' and n, we get C(z') > C(u) — 
2w max = w* — 2w max . 

By Lemma [TT| matching z" satisfies t T z" < t T v and w T z" > w T v — 4w max . We observe that z" 
satisfies the budget constraints since 

fz" < £ T v = £ T {( ai + a 2 )z' + a 3 x 3 ) < £ T ({ ai + a 2 )u + 03^3) = £ T x* < L. 

Furthermore, 

w T z" > w T v - 4w max = w T ((ai + a 2 )z + a 3 x 3 ) - 4w max 
> w t ((qi + a 2 )u + a 3 x 3 ) - 6w max = w* - 6w max . 

Let OPT' be an optimum solution to the reduced problem. Of course, w* > w(OPT'). Further- 
more, the weight of the guessed edges Eh is at least Q/ew max . Since w(OPT) = w(Eh)+w(OPT'), 
we can conclude that the solution returned by the algorithm has weight at least w(£j^)(l — e) + 
w(OPT') > (1 - e)w(OPT). □ 

5 Conclusions 

A first obvious open problem is finding a PTAS for the fe-BUDGETED MATCHING problem for any 
k = 0(1). It is interesting to notice that most parts of the approach presented in Section [J] can easily 
be generalized to an arbitrary constant number of budget constraints. More precisely, the only part 
that is tailored to two budgets is Lemma [71 which is only valid for curves in two dimensions. We 
believe that for every constant k € N there is r(k) € N such that the following generalized version 
of Lemma [7] holds. 

Conjecture 12. Let f : [0,r] — > M. k be a curve, and let \x € [0,r]. Then there are r(k) disjoint 
intervals [a%, bi], . . . , [a r(fe ), 6 r(fc) ] C [0, r] such that f '(0) + X)i=i (/ V>i) ~ f ( a i)) = + i 1 ~ v)f( T )- 
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In particular, for k = 2, the pair (a,t) given by Lemma [7] defines either one interval [a, a + t] or two 
intervals [a, r] and [0, a + t — r] satisfying the claim above. If Conjecture [12] holds, then a PTAS 
for the fc-budgeted case can be obtained by following a procedure analogous to the one that we 
presented here. 

Another interesting direction for further research, is to check whether a similar technique can be 
applied to ^-budgeted matroid intersection independent set. A crucial property that 
we exploited in this paper is the well-known fact that a matching can be transformed to another 
matching in the same graph by exchanging edges on alternating paths and cycles. Similar exchange 
properties are known for the intersection of matroids. However, they do not seem to allow for an 
easy adaption of the presented algorithm. 

Another set of questions involves the problems considered here but with one budget. Is there a 
fully-polynomial PTAS (FPTAS) for 1-budgeted spanning tree and 1-budgeted matching? 
We remark that, as noted in [6], an FPTAS for the second problem would imply a deterministic 
algorithm for EXACT MATCHING with polynomial weights, which is a long-standing open problem. 
In the case of 1-budgeted perfect matching, a PTAS is not known. 



Acknowledgements. The authors wish to thank Friedrich Eisenbrand for supporting their visit 
to EPFL and for helpful discussions. 
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Figure 1 The full line indicates /, and the dashed lines f ai and f a2 , respectively. 




Appendix 

Proof of Theorem We show that deciding feasibility of the considered problems is iVP-complete. 
It is sufficient to prove the claim for k = 2. Consider first 2-budgeted spanning tree: the claim 
for A;-BUDGETED MATROID BASIS and, consequently, for A;-BUDGETED MATROID INTERSECTION 
BASIS trivially follows. Let T >+ denote the problem, and V^ 1 its variant with arbitrary (i.e., positive 
and/or negative) lengths. Of course, includes V + as a special case. To see the opposite 
reduction, observe that a spanning tree contains exactly n — 1 edges. Hence, by adding a sufficiently 
large value M to all the lengths, and adding (n — 1)M to the budgets, one obtains an equivalent 
problem with non- negative lengths. It is easy to see that includes as a special case the problem 
V = of determining, for a given length function £'{■) and target L', whether there exists a spanning 
tree T of length £'(S) = L': a reduction is obtained by setting ^i(-) = — = £'(■) and Li = 
— L2 = L'. Hence it is sufficient to show that V = is NP-complete. We do that via the following 
reduction from partition: given a\, a%, . . . , a q £ Q and a target A £ <Q, determine whether there 
exists a subset of ctj's of total value A. Consider graph G q , consisting of q cycles C\, C2, ■ ■ ■ , C q , 
with C{ = (flj, bi, Ci, di) and q = aj+i for i = 1, 2, . . . , q — 1. Let I'iaibi) = ati, i = 1, 2, . . . , k, and 
set to zero all the other lengths. The target is L' = A. Trivially, for each spanning tree T and each 
cycle C{, the length of T n Cj is either or «j. Hence, the answer to the input partition problem 
is yes if and only if the same holds for the associated instance of V = . 

Consider now 2-budgeted perfect matching. Since each perfect matching contains exactly n/2 
edges, with the same argument and notation as above it is sufficient to prove the ./VP-completeness 
of the problem V = of determining, for a given length function £'(■) and target L', whether there 
exists a perfect matching M of length £'(M) = L'. We use a similar reduction from partition as 
above. The graph is again given by the cycles C\, . . . ,C q . However, this time each cycles forms 
its own connected component. Furthermore, the lengths are given by £'{aibi) = £'{cidi) = CKj, 
i = 1, 2, . . . , k, all the other lengths are zero, and L' = 2A. It is easy to see that, for each perfect 
matching M and each cycle Q, the length of T n is either or 2a j. The claim follows. 
Eventually consider 2-budgeted shortest path. We restrict our attention to the graph G q as 
used for the spanning tree reduction, and let (s,t) = (a±,c q ). Since any s-t path in this graph uses 
exactly 2q edges, we have by the usual argument that it is sufficient to show the ./VP-completeness 
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of the problem V = of determining, for a given length function £'(■) and target L', whether there 
exists an s-t path P of length £'{P) = L'. The claim follows by essentially the same reduction as 
in the spanning tree case. 

□ 

Proof of Lemma^ Without loss of generality, we can assume that / = (/i, $2) satisfies. 

(i) /(o) = (0,0) and /(r) = (r,0) for r > 0. 

(ii) / 2 (t)>0Vie[0,r], 

(iii) / is not self-intersecting (i.e., / is an injection). 

Let / = (/i,/2). Since the statement of the lemma is independent of changes in the coordinate 
system (and the claim is trivial for /(0) = /(t)), we can assume that /(0) = (0, 0) and /(r) = (r, 0) 
for some r > 0. Hence, (i) holds. The Gasoline Lemma [6] states that there is a € [0, r] such that 
f 2 (t) > Vi G [0,t]. More precisely, this condition is satisfied by choosing a G argmin{/2(i) | t G 
[0,r]}. One can easily observe that if the lemma is true for f a then it also holds for /. Hence, 
we can assume, by replacing / by f a , that Property (ii) holds. Property (iii) can be enforced by 
removing loops. 

Furthermore, we assume [i G (0, 1), otherwise the claim is trivially true. Let v = («i,0) = nf(0) + 
(1 — p)f{r) = (1 — p)/(r) and g : [0, r] — > M 2 be the translation of / by v. 

g{t) = f{t)+v VtG[0,T]. 

Let u = max{f 2 (t) | t G [0,r]} and p = min{i G [0,r] | /2(i) = u}, i.e., /2 attains its maximum 
the first time at p. Let /' : [p, t] — > M 2 with /'(i) = /(t) be the subcurve of / over the interval 
[p, t], and let </ : [0,p] — >■ M 2 with ^'(t) = g(t) be the subcurve of g over the interval [0,p]. In the 
following we show that /' and g' intersect. Consider the endpoints of /' and g' . The endpoints 
/'(r) = (r, 0) and ^'(0) = (v±, 0) both lie on the x-axis, and since v± = (1 — /u)r < r, /'(r) lies to the 
right of ^'(O). Similarly, the other two endpoints f'(p) = (fi(p),u) and g'(p) = (flip) + v\,u) have 
the same y-components, where this time f'(p) is to the left of g'(p). Since the second component of 
both curves lies between and u, one can easily deduce that they have to cross. In more detail, one 
way to show this is to consider the polygonal curve h joining, in the given order, f'(p) = (fi(p),u), 
(— M, u), (— M, —1), (r, —1) and (r, 0) = f'(r), where M > is a large value such that h does not 
intersect g' . The concatenation of h with f forms a closed, not self-intersecting curve, which - by 
Jordan's Curve Theorem [22J - divides M? into two regions, a bounded one and an unbounded one. 
Furthermore, g' has one endpoint in one region and the other endpoint in the other region. Again 
by Jordan's curve theorem, we have that g' and /' intersect (since we have by construction that g' 
does not intersect h). Hence, there exists t\ G [p, r] and t% G [0,p] such that f(t\) = gfo). Since 
ftp) 9ip)i we have ti <t\. The claim is satisfied by 

f t2 (h - t 2 ) = f(tx -t 2 + t 2 ) - f(t 2 ) = g(t 2 ) - f{t 2 ) = v. 

□ 
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